Quantum Measurements of Coupled Two-Level Systems 
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We propose an approach to measuring nonresonant coupled systems, which gives a parametrically 
smaller error than the conventional fast projective measurements. The approach takes into account 
that, due to the coupling, excitations are not entirely localized on individual systems. It combines 
high spectral selectivity of the detector with temporal resolution and uses the ideas of the quantum 
diffusion theory. The results bear on quantum computing with perpetually coupled qubits. 
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I. INTRODUCTION 

The understanding of quantum measurements has sig- 
nificantly advanced in recent years, in part due to the 
fast development of quantum information theory [l|, Q • 
Measurements constitute a necessary part of the oper- 
ation of a quantum computer. In the context of quan- 
tum computing, it is often implied that measurements 
are performed on individual two-state systems, qubits 
and that during measurements qubits are isolated from 
each other. However, in many proposed implementations 
of quantum computers the qubit-qubit coupling may not 
be completely turned off. The interest in measuring cou- 
pled systems is by no means limited to quantum com- 
puting; however, qubits provide a convenient language 
for formulating the problem. 

In a system of coupled qubits (coupled quantum sys- 
tems) excitations are not entirely localized on individ- 
ual qubits even where the qubits have different energies. 
Therefore a one-qubit measurement may miss the exci- 
tation mostly localized on the qubit which is being mea- 
sured. A measurement may also give a false-positive re- 
sult: the detected excitation may be mostly localized on 
another qubit but have a tail on the measured qubit. In 
the context of quantum computing, this is a significant 
complication, since the overall error accumulates with the 
number of qubits. 

The most simple type of measurements which was es- 
sentially implied above is a fast von Neumann-type pro- 
jective measurement Q. Alternatively one can use con- 
tinuous measurements, in which the signal from the qubit 
is accumulated over time ■ They are often implemented 
as quantum non-demolition measurements (QNDMs) in 
which the quantity to be measured (like population of 
the excited state) is preserved while a conjugate quan- 
tity (like phase) is made uncertain. As we will see, the 
standard QNDMs do not solve the precision problem for 
interacting qubits. 

The goal of the present paper is to find a way of mea- 
suring nonresonantly coupled qubits that gives a para- 
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metrically smaller error than standard fast or continuous 
one-qubit measurements. The idea is to combine tem- 
poral and spectral selectivity so as to take advantage of 
different time and energy scales in the system. The pro- 
posed measurement is continuous. However, it is not of 
a QNDM-type, both the qubit population and the phase 
of the qubit wave function arc changed. 

We will consider the situation where excitations are 
strongly (albeit incompletely) localized on individual 
qubits, at least for the time that significantly exceeds 
the duration of a measurement. This can be accom- 
plished by appropriately tuning the qubit energies e„. 
For one excitation, localization of stationary states is 
well-understood since Anderson's work Q on disordered 
systems where e n are random. Anderson localization re- 
quires that the bandwidth h of the energies e n be much 
larger than the typical nearest-neighbor hopping integral 
J. One-excitation localization becomes stronger for the 
same hj J, i.e., the localization length becomes smaller if 
£ n are tuned in a regular way so as to suppress resonant 
excitation transitions. The problem of localization of 
multiple excitations is far more complicated because the 
number of states exponentially increases with the number 
of excitations. However, at least for a one-dimensional 
qubit system, by tuning e n one can still obtain strong lo- 
calization of all excitations for a time cx J _1 (/i/ J) 5 @,@]; 
here and below we set h = 1 . 

We propose to detect an excitation by resonantly cou- 
pling the measured qubit (MQ) to a detecting two-level 
system (DS). If the MQ was initially excited and the DS 
was in the ground state, the excitation can move to the 
DS. There its energy will be transferred to the reservoir, 
and the change of the state of the reservoir will be di- 
rectly detected. For example, the DS can emit a photon 
that will be registered by a photodetector. The typical 
rate of photon emission T should largely exceed the rate 
of resonant (but incoherent) excitation hopping from the 
MQ to the DS, so that the probability for the excitation 
to go back onto the MQ be small. The rate T should 
also largely exceed the interaction-induced shift of the 
MQ energy levels. At the same time, T should be small 
compared to the bandwidth of site energies h. Then the 
qubits adjacent to the MQ are not in resonance with the 
DS, and the rate of excitation transfer from these qubits 
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to the DS is small. 

If the above conditions are met, there should be a time 
interval within which an excitation localized mostly on 
the MQ will be detected with large probability, whereas 
excitations localized mostly on neighboring qubits will 
have a very small probability to trigger a detection signal. 
As we show, the associated errors are much smaller than 
in a projective measurement. 

State measurements for coupled qubits have been per- 
formed with Josephson-junction based systems, where 
there were studied oscillations of excitations between 
the qubits 0, @|. The oscillations occurred where the 
qubits were tuned in resonance with each other. The 
effect of the interaction on measurements in the case 
of detuned qubits, which is of interest for the present 
work, was not analyzed. The approach proposed here, 
which combines spectral and temporal selectivity to en- 
able high resolution is different not only from the stan- 
dard continuous quantum measurements 0] but also from 
other types of time-dependent quantum measurements, 

cf. Rcfs. Mosul 

In Sec. II we describe the system of two qubits and a 
resonant inelastic-scattering based DS. We identify the 
range of the relaxation rate of the DS and the qubit pa- 
rameters where the measurement is most efficient. In 
Sec. Ill the description of the time evolution of the sys- 
tem is given and it is shown that the decay rates of 
the stationary states of the qubits are strongly different. 
The detailed theory for one- and two-qubit excitations is 
given in Appendices [A] and [Bj respectively. Section IV 
describes how the initial states of the system can be effi- 
ciently discriminated, and the analytical results arc com- 
pared with a numerical solution. Section V contains con- 
cluding remarks, including an extension of the results to 
a multi-qubit system. 



II. THE MODEL 

A. A two-qubit system 

We will concentrate on a quantum measurement of two 
coupled two- level systems (two spin- 1/2 particles or two 
qubits) and then discuss how the results extend to a 
multi-qubit system. The system is sketched in Fig. [1] 
We will assume that the qubit excitation energies (the 
spin Zeeman energies) £i ; 2 largely exceed both the inter- 
action energy J and the energy difference |tL>2i|, where 
^21 = £2 — £i; for concretcness, we assume that u>2i > 0. 
Via the Jordan- Wigner transformation the system can be 
mapped onto two spinlcss fcrmions with the Hamiltonian 



fermions present. A fermion on site n corresponds to 
the nth qubit being excited. The parameter J in this 
language is the hopping integral, whereas J A describes 
the interaction energy of the excitations. 

We will assume that J <C u>2i- In this case the station- 
ary single-particle states of the system l^i 2 ) are strongly 
localized on sites 1 or 2, 

|V>< rt) >«C(|n>-/i(-l)"|3-n>) (n = l,2), (2) 

where fi = 2fe 1 /J « -J/2w 2 i and C = (1 + ^ 2 )" 1/2 
(&l = [ujix — V '^21 + J 2 ]/^ i s the shift of the energy 
level £1 due to excitation hopping). 




2nd Qubit 



1st Qubit Detecting System 
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FIG. 1: (Color online) The measurement scheme. Qubits 1 
and 2 are not in resonance, but are perpetually coupled, with 
the coupling constant J <C ^21, where W21 = £2 — £1 is the 
level detuning, 1021 <C £1,2- The detecting system is reso- 
nantly coupled to qubit 1. When an excitation is transferred 
to the DS, the DS makes a transition to the ground state (for 
example, with photon emission) which is directly registered. 
The transition rate F exceeds J, Jd , but is small compared to 
W21 to provide spectral selectivity. 



Measurements are done by attaching a detector to a 
qubit, i.e., to the physical system represented by the 
qubit. For concretcness, we assume that the measured 
qubit (MQ) is qubit 1. It is seen from Eq. @ that, if 
the measurement is fast projective and the system is in 
state {ip^}, the occupation of this state will be detected 
with an error fi 2 . With probability fj? the detector will 
"click" also if the system is in state IV'^*^)- The same 
error occurs, evidently, in a simple continuous measure- 
ment that does not involve energy transfer, like measure- 
ments with quantum point contacts or tunnel junctions 
[HI, [l^. Moreover, the decoherence of the qubit brought 
about by such a measurement, even where its rate is small 
compared to w 2 i, will ultimately lead to spreading of the 
excitation over both states ) , which leads to an extra 
limitation on the measurement precision. 



B. Resonant inelastic-scattering detector 



Here, the subscript n = 1,2 enumerates the coupled 
qubits, and a n ,ajj are the fermion annihilation and cre- 
ation operators. The ground state corresponds to no 



The state of the system can be determined with a 
higher precision using a detector that involves inelastic 
transitions. A simple model is provided by a two-level 
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DS which is resonant with qubit 1. sec Fig. [T] The mea- 
surement is the registration of a transition of this system 
from its excited to the ground state; for example, it can 
be detection of a photon emitted in the transition. The 
Hamiltonian of the qubit-DS system is 

H = H s + £Da D a D + -J D (a{a D + a D a>i), (3) 

where Ed is the energy of the excited state of the DS 
and Jd characterizes the coupling of the DS to qubit 1, 
J d < w 2 i. 

The qubit-DS dynamics can be conveniently analyzed 
by changing to the rotating frame with a unitary trans- 
formation U(t) = exp[— ie±t J2a=i 2 d a a a a\- We assume 
that relaxation of the DS is due to coupling to a bosonic 
bath (photons). Provided this coupling is weak, the 
qubit-DS dynamics in slow time (compared to e^ 1 ) is 
described by a Markov equation 



p = i[p, H]-T [a D a D p - 2a D pa D + pa D a D 



(4) 



Here, H = H — e\ J2 a =i 2 d a a a a, an d T = T(ed) is the 
DS decay rate. Equation ([4]) applies provided T <C £1; 
in addition, the dispersion of T has been disregarded, 
\uj2i(dT / de d)\ *C T. The rcnormalization of Ed due to 
the coupling to the bath is assumed to be incorporated 
into ££>, and the bath temperature T <C eo/ks, so that 
there are no spontaneous transitions of the DS from the 
ground to the excited state, 

Equation (fj| should be solved with the initial condition 
that for t = the DS is in the ground state whereas the 
qubits are in a state to be measured. As a result of the 
excitation transfer from the qubits, the DS can be excited 
and then it will make a transition to the ground state. 
The directly measured quantity is the probability R(t) 
that such a transition has occurred by time t, 



R(t) = 2T dtTr p{t)a D a D 



(5) 



It is clear, in particular, that if one of the qubits is ex- 
cited, the excitation will be ultimately fully transferred 
to the DS and then further transferred to the photon 
bath, so that R(t) — ► 1 for t — > 00. If on the other hand 
both qubits are in the ground state, then R(t) = 0. 
Qubit measurements can be efficiently done for 



w 2 i > T > J, Jd, \e d - ei\ 



(G) 



In this case, there are no oscillations of excitations be- 
tween qubit 1 and the DS. The energy detuning Ed — £\ 
plays no role, and without loss of generality we can set 
Ed = £i- We will start the analysis with the case where 
there is no more than one excitation on the qubits for 
t = 0. Then one can replace H with a single-excitation 
Hamiltonian, 



H 



w 2 ia 2 a 2 



2 ft 2 + 7: (ja\a,2 + JDa\ a D + Ha 



As we show, for a DS that satisfies conditions (J6)) time 
evolution of R(t) is characterized by two strongly differ- 
ent scales, which makes it possible to differentiate be- 
tween the qubit states. 



III. TIME EVOLUTION OF THE DENSITY 
MATRIX 

A. One excitation 

The continuous measurement ([5]) requires finding ex- 
pectation values (at(t)a a (t)) = Tra^a a p(t) = p a p(t), 
where a,/3 run through the subscripts 1,2,_D. The ma- 
trix elements p a p = p* ia satisfy a system of nine lin- 
ear equations that follow from the operator equations 
dU), (J7J). This system of equations is closed, the ma- 
trix elements p a p do not mix with the expectation values 
(o a (t)), (a>a(t))- In the case where there is only one exci- 
tation (one fermion) they also do not mix with expecta- 
tion values (a' r a (t)ap(t)a.y(t)as(t)) , which are all equal to 
zero. 

The solution of the equations for Pap{t) and the anal- 
ysis of the signal R(t) are simplified in the range ©. 
The relaxation rate of the matrix elements po a that in- 
volve the DS is T. For T ^> Jd, J this rate is faster 
than other relaxation rates, as explained below (see Ap- 
pendix A for details), and therefore over time the 
matrix elements pDa reach their quasi-stationary values. 
Relaxation of pn, on the other hand, is determined by 
the excitation transfer from site 1 to the DS. This trans- 
fer is similar to quantum diffusion of defects in solids 
within narrow bands for weak coupling to phonons fl3| 
or between closely spaced discrete non-stationary states 
of weakly coupled reorienting defects It is charac- 

terized by rate 



Wi = j 2 d /2y, Wi < r. 



(8) 



Equation JS} can be readily understood in terms of the 
Fermi golden rule: this is a transition rate from qubit 1 to 
the DS induced by the interaction oc Jd , with T being the 
characteristic bandwidth of the final states and T _1 being 
the density of states at the band center, respectively. The 
rate W\ describes relaxation of an excitation localized 

initially in state lip^)- 

If the excitation is localized mostly on qubit 2 and 
has energy ks aj 2 i, its decay rate becomes much smaller. 
The decay can be thought of as due to the interaction of 
qubit 2 with the DS, which is mediated (nonresonantly) 
by qubit 1. Therefore the effective interaction energy 
is w JJd/2uj21- The effective density of final states is 
determined by the tail of the density of states of the DS 
at frequency LO21 S> T. For exponential in time relaxation 
described by Eq. (0| the density of states is Lorentzian, 
and on the tail at frequency u>2i it is ~ r/wfi- Therefore 
the expected decay rate W2 is 



(7) 



W 2 = J a J d T/8uj^ W 2 < Wi. 



(9) 
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It follows from the result of Appendix [5] that, if for 
t = the two-qubit system is in a stationary state l^n^} 
given by Eq. flSJ) (ra = 1,2), the probability R„(t) to 
receive a signal by time t is 



R n (t) = 1 - exp(-W„t) (n=l,2). 



(10) 



The strong difference between and T justifies the 

assumption that po a reaches a quasistationary value be- 
fore the populations pn, P22 change. The full adiabatic 
theory of qubit relaxation is described in Appendix A. 
We note that the analysis takes into account corrections 
~ (J/uj2i) 2 , which are important, since they determine 
the difference between a fast projective measurement and 
the continuous measurement discussed here. 



B. Two excitations 

The above analysis can be readily extended to the case 
where both qubits 1 and 2 are initially in the excited 
state, i.e., initially there are two spinless fcrmions on 
sites 1 and 2. We will assume that the interaction be- 
tween the excitations (fcrmions) J A is not strong, so 
that |JA| <C r. The qualitative picture of the dynamics 
of the system then is simple. First, over time ~ T _1 there 
is established a quasi-stationary "quantum diffusion cur- 
rent" from qubit 1 to the DS, which is determined by the 
transition rate W\ and is equal to Wipn{i). It drains 

stationary state |'0 1 st ' ) ) over time ~ W-f x . The presence 

of an excitation (fermion) in state |V*2 ) on ly weakly af- 
fects this current because of the large energy difference of 
the states u>2\- After state l^i ) i s emptied, further evo- 
lution corresponds to the single-excitation decay of state 
|^(st)^ Therefore the overall signal should be 



R (2e \t) = R 1 {t)+R 2 {t) = 2-e 



-Wit 



,-Wat 



(11) 



A detailed derivation of this expression is given in Ap- 
pendix [B] It is clear from the above analysis that the 
most interesting problem is to distinguish which of the 
stationary states |V ; i2 } is initially occupied; the situation 
where both of them are occupied is simpler. 



IV. RESOLVING ONE-EXCITATION STATES 
FROM TIME-DEPENDENT MEASUREMENTS 

The rates W\ and W2 of signal accumulation for differ- 
ent initially occupied one-excitation states are paramctri- 
cally different. This enables efficient distinction between 
the states using a resonant DS. The results for time evolu- 
tion of the signals R(t) obtained by a numerical solution 
of the system of equations for the matrix elements p a p 
are shown in Figs. [2] and [3J The figures refer, respec- 
tively, to the cases where the system of qubits is initially 
in one of the eigenstates 1^12^) and in a mixed state. It is 



seen from the figures that, even where the parameter ra- 
tios T/LJ21 and J/T, J/j/r are not particularly small, the 
results are well described by the asymptotic expressions 

(nsi. 




FIG. 2: (Color online) Time dependence of the probability to 
detect a signal for the two-qubit system being initially in the 
stationary states |'i/>i St ' 1 ) (a) and IV4 ) ( D )- Note the differ- 
ence of time scales. The parameters are UJ21/T = 4, J/T = 
Jd/T = 1/2. The solid and dashed lines show, respectively, 
the numerical solution of the master equation and the asymp- 
totic expressions (|10|l for the limiting case 0J21 > T > J,Jd- 
The insets show | log[l — -Ri (*)] | for small time, where the dif- 
ference between the numerical and asymptotic expressions is 
most pronounced. 



Expressions (|10|) are obtained in the adiabatic approxi- 
mation and do not describe how the matrix elements p a p 
reach their adiabatic values. In particular, they do not 
describe the evolution of R for t < Since the DS 

is initially in the ground state, it follows from Eqs. (j4|) 
and ([5]) that R(t) oc t 2 for t — > in contrast to R cx t as 
predicted by Eq. iflO]). 

Breaking of the adiabaticity for short times explains 
the shifts of the asymptotic curves (jTTJJ) with respect to 
the numerically calculated curves in Fig. [2] A simple 
estimate for the shift can be obtained for detection of 
state lip^} ■ Here, adiabaticity is established over time 
t ~ r _1 ; this is the relaxation time of the matrix ele- 
ments p a D, PDa that directly involve the DS. Then from 
Eq. (fTOf the shift should be ~ Wi/T, which agrees with 

Fig. [21a). For the initially occupied state {ip^) * ne ex " 
ponential decay of i?2, Eq. (flT))) . is obtained in the doublc- 
adiabatic approximation, which exploits the interrelation 
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between the relaxation rates T 3> W\ 3> W?. The adia- tive measurement on qubit 1 gives 



batic regime is formed over time ~ W{~ , and the shift of 
the numerical curve with respect to the asymptotic one 
is ~ W2/W1, which agrees with Fig. [21(b). 




FIG. 3: (Color online) Time dependence of the probability 
to detect a signal for the two-qubit system being initially in 
a superposition of stationary states \tp(0)} = cos 9\ip[ st ' 1 ) + 
sin 8\i))^ ') . Curves 1 and 2 present the numerical solution of 
the master equation and refer to 9 = 7r/4 and n/3, respec- 
tively. The parameters of the qubits and the DS are the same 
as in Fig. [2] Inset: the solid line shows R(t) for 9 — ±7r/4 
(Pi = 1/2) in the optimal time range (|14| l. the dashed lines 
show the results of fast projective measurements on qubit 1, 
|<lhA(O)}| 2 ,for0 = ±7r/4. 



The vast difference of the time scales over which the 
signal accumulates depending on the initially occupied 
state allows one to discriminate between the states with 
high precision. If the initial state of the system | -0(0)) 

is a superposition of stationary states and IV 7 ^); 

for Tt ^> 1 the signal R(t) is the appropriately weighted 
superposition of the signals Ri(t) and R^it), 

R(t) w Pi [1 - exp(-Wit)] + P 2 [1 - exp(-W a *)] , 

p„ = |<v(oM st) >f- (12) 

Here, Pi^ are the initial populations of the stationary 
states. We note that Pi + P2 ^ 1, in the general case. 
The initial state of the two-qubit system can be a super- 
position of the ground and one- or two-excitation states. 

Over time ~ W^ 1 the function R(t) approaches the 
population Pi of state {tp^). The further change of R(t) 
occurs over a much longer time ~ W^ ■ This is seen in 
Fig. [3] As explained in Appendix [XJ a contribution to 
R(t) from the fast-oscillating [as exp(±iu>2it)] terms in 
the density matrix is small compared to {J/0J21) 2 : for 
Tt ^> 1 the corresponding correction to Eq. (fT2")) is < 
J J"jjjAu)\ 1 Repi2(0) <C J 2 /4w2i (we assume Jd ~ J)- 

From Eq. (|12|) . in a broad time interval the error in 

the measured population of state |'0 1 st '') is smaller than 
in a fast projective measurement. If the initial state is 



|<1|V(0)>| 



J 2 



P 



1/2 



P 



2w 2 i 



(13) 



This differs from Pi by ~ J/u>2i f° r Pi ~ P2; m the case 
of strongly different populations where P1/P2 or P2/P1 
is i5 J 2 l u >2\ the difference becomes ~ J 2 /cjfi- The pro- 
posed measurement gives a much smaller (in fact, a para- 
metrically smaller) error. We have 



\R{t)-Pi\ « 



J 



4^ 2 2 i 



for 



3 -Wit 



,W 2 t < 



J 2 

4^21 ' 



(14) 



V>(0)) = P x /2 \il>[ st) ) + P 2 1/Z exp(^)|^r ; ). a fast projec 



,1/2 



(st) N 



In the explicit form, the time interval for high-accuracy 
measurement of state |i/ , i St ^) is determined by condition 
21n(2w 2 i/J) < Wit < (w 2 i/r) 2 . This condition is easy 
to satisfy in the range ([6]). We emphasize that, as seen 
from Eq. (fl~3"]) . in most cases the error of a fast projective 
measurement is parametrically larger than J 2 /ui^i- The 
difference between the proposed approach and a fast pro- 
jective measurement is illustrated in the inset in Fig. [3] 
Advantageous features of the proposed approach can 
be seen from a comparison with a similar but a seem- 
ingly simpler scheme in which one directly turns on the 
coupling of qubit 1 to a thermal reservoir and detects 
the emitted excitation. This scheme is analogous to the 
scheme implemented in Josephson phase qubits @, Q 
where decay of one of the qubits was effectively "turned 
on" by reducing the appropriate tunnel barrier. In such 
a scheme the decay rate of the initial stationary state 
localized mostly on qubit 1 is determined by the intro- 
duced relaxation rate of this qubit T. The decay rate 
of the stationary state localized mostly on qubit 2 can 
be estimated, following the arguments of Sec. IIII1 as 
~ J 2 r/w| 1; provided w 2 i > T,J. Therefore, if the 
system is initially in state \4>2), over the lifetime T _1 
of state |f/'i St ' ) ) the detector will "click" with probability 
~ (J/tj2i) 2 rt > (J/UO21) 2 ■ As a result, there is no gain 
compared to a fast projective measurement. 



V. CONCLUSIONS 

In this paper we have demonstrated that a state of two 
coupled two- level systems (qubits), which are detuned 
from each other, can be determined with an accuracy 
much higher than that of a fast projective measurement. 
This is accomplished by combining frequency and tem- 
poral selectivity. The detector is tuned in resonance with 
the qubit to which it is directly coupled. When the detec- 
tor relaxation rate Y is higher than the coupling J, Jd but 
less than the energy difference between the qubits UI21, 
the overall relaxation of the system is characterized by 
two strongly different time scales. One relaxation rate, 
W\ = Jp/2T, determines the decay of the excited sta- 
tionary state localized mostly on the resonant qubit. The 
decay rate of the excited stationary state localized mostly 
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on the non-resonant qubit W2 = (J Jd) 2 ^/ 81^21 is para- 
metrically smaller. This makes it possible to discriminate 
between the states using time-dependent measurements. 

The measurement requires turning on the interaction 
between the measured qubit and the DS. Alternatively, 
and often more conveniently, one can tune the DS in res- 
onance with the qubit. The required precision is deter- 
mined by the decay rate of the DS T. Moreover, one can 
slowly sweep the DS energy level through the qubit en- 
ergy, so that the energies stay in resonance (to accuracy 
~ T) for a time > Wf 1 . 

The results can be generalized to the case of a many- 
qubit system. Of particular interest is a qubit chain with 
Hamiltonian 

H = ^2 £ nala„ + - J ^(a^a„+i + a ] n+l a n ) 

n n 

+JA^a, t l a r t l+1 a„+ia„, (15) 

n 

which is an immediate extension of the two-qubit Hamil- 
tonian (fT]). As mentioned in Sec. I, for appropriately 
tuned site energies e n all many-particle excitations in 
such a chain remain localized for a long time, which scales 
as a high power of h/J (h is the typical bandwidth of 
energies e n ) d, 0]. However, excitations are not fully 
localized on individual sites, in (quasi)stationary states 
the tails of their wave functions on neighboring sites are 
~ J/h. 

Following the proposed method, to determine whether 
there is an excitation localized in a quasistationary state 
centered at a given qubit one should couple this qubit 
resonantly to a DS. The excitation will be detected over 
time ~ . Excitations localized mostly on neighbor- 
ing qubits will not affect the measurement as long as 
the duration of the measurement is small compared to 
W% • The interaction between excitations will not affect 
this result for appropriately chosen decay rate of the DS. 
The reduction of the measurement error is an important 
prerequisite for scalable quantum computing with per- 
petually coupled qubits. 



APPENDIX A: THE ADIABATIC 
APPROXIMATION 

In the presence of one excitation on the qubits, the 
dynamics of the two-qubit system coupled to the DS 
is described by Eqs. (H}, ([7]), which can be written as 
nine linear first-order equations for the matrix elements 
Pap(t) = (ap(tya a (t)). Formally one can solve these 
equations by finding the corresponding eigenvalues and 
eigenfunctions. For a strong inequality between the tran- 
sition frequency W21 , on the one hand and, on the other 
hand, the decay rate T and the hopping integrals J, Jd 
©; the eigenvalues can be separated into those corre- 
sponding to fast weakly damped oscillations and to slow 
evolution. There are 4 "fast" eigenvalues with imaginary 



part close to ±W2i ■ In particular, as seen from Eq. ([4]), 
to zeroth order in J, Jd we have pn^poi °c exp(iW2ii)- 
In turn, where T J, the "slow" eigenvalues can be 
separated into three eigenvalues with real parts cx — T. 
In particular, to zeroth order in J, Jd we have pdd 
exp(— 2r<) and poi = p\ D oc cxp(— Tt) . Then there arc 
two more "slow" eigenvalues, which as we will show are 
given by Eqs. © and ©. 

We start with the effect of slow motion on the mea- 
sured signal R(t). The analysis of slow dynamics can be 
done in the adiabatic approximation. The relaxation rate 
of the matrix elements pr>a {pt = 1, 2, D) is oc T. There- 
fore over time t ^> these matrix elements reach quasi- 
stationary values. These values adiabatically follow the 
slowly evolving matrix elements p nm , where Roman sub- 
scripts n, m enumerate the excited states of the qubits. 
In particular, by noting that (d/dt) ^ Q p aa = —2TpoD 
we obtain 

Pdd » -(pu +p 22 )/2r. (Al) 



1. Evolution of the one-excitation stationary state 
resonant with the DS 

We assume first that the system is in the stationary 
state \ip^), which has energy close to the DS energy. In 
this case, for t = we have pu P22 ~ (J/^2i) 2 pu 
and all off-diagonal matrix elements p a p with a ^ (3 are 
small compared to pu. As we will see, this hierarchy per- 
sists for t > as well. For t ^> to leading order in 
J, J D we have p D1 w i(J D /2T)(p DD - p n ) + i( J/2T)p D2 
and pdd ~ — (JD/2r)Impui. Substituting pui into the 
equations for pu and eliminating pjjjj we obtain, to lead- 
ing order, 

pii«-Wipii+Aii, (A2) 

All = 7}iJ(pi2 - P2l) + 4^T-RcpD2, 

where the relaxation rate W\ is given by Eq. |5]). In 
agreement with the adiabaticity assumption, we have 
Wi <C r. The term An docs not affect the evolution 
of pu if the system is initially in state \ip^}, but it is 
important for the analysis of evolution from state {ip^)- 
To find time evolution of P22 one has to use equations 

P12 = i^2iPi2 + -^iJ(Pn - P22) - ^JdPd2, (A3) 

P22 = --jiJ{pi2 - P21) 

that immediately follow from operator equation . Us- 
ing the similar equation for f>£>2 and the expression for 
Pdi given above one can show that, for pn ^ P22, 
PD2 ~ i(JJD/4u>2iT)pii. This expression should be 
substituted into Eq. (|A3|) for pi 2 . One than obtain 
Impi2 ~ (J/4w|i)pn. Then from the second equation 
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[3]p22 ~ (•//2w2i) 2 pii- Therefore, if the system is ini- 
tially in the stationary state lip^}, the relation between 
P22{t) and p\i(t) remains unchanged, to leading order in 
J. 

From Eqs. ([AT]) . (|A2j) we see that, if the two- 

qubit system is initially in the stationary state {ip^}, 
the probability to have detected this state by time £ is 
Ri(t) = cxp (—Wit). We use this expression in Eq. (fTU)) . 



2. Evolution of the nonresonant one-excitation 
stationary state 

If the system is initially in state we have p\\ <C 

P2i- Then it is important to keep track of higher-order 
corrections in J/w 2 i in the equations for p a p. In par- 
ticular, in the adiabatic approximation one has to write 
PD2(t) as 



PD*{t) » [2(w2i+*T)]- 



(A4) 



The decay rate of Pi2{t) can be obtained by substi- 
tuting Eq. (|A4[) into Eq. (|A3[) for pi2, which gives for 
the decay rate an expression (Jd /2w2i) 2 r '. This rate is 
small compared to the decay rate T of pD2, justifying the 
adiabatic approximation used in Eq. (|A4[) . At the same 
time, it is large compared to the decay rate of p22- To 
describe decay of P22 , we substitute the adiabatic solu- 
tion of Eqs. (|A3[) for Impi 2 [with account taken of (|A4[) ] 
into Eq. (|A3[) for p 2 2- This gives 



P22 ~ -W2P22, 



(A5) 



with the decay rate W2 of the form of Eq. © . 

Noting that Repi2 « (J/2w2i)p22) we obtain from 
Eqs. (|A"2j) . (TM| that, for time £ > Wf 1 , we have 
pn(t) « (J/2w 2 i) 2 /9 2 2(i). Then from Eqs. (0 (pO}, and 
(|A5|) we obtain the explicit expression (fT0|) for the prob- 
ability to detect state |V4 ) m time £. 

If initially the system is in a superposition of the 
states Itp^}: along with slowly varying (on time u^i) 
solutions there emerge terms in the density matrix that 
contain fast oscillating factors oc exp(±z<x>2i£)- One 
can show that, for £ > their decay is con- 

trolled by the decay of the fast-oscillating component 

P12 W = [P2i(0]*i which is characterized by the fac- 
tor exp(— W\t/2). The corresponding fast-oscillating 
term in the population of the excited state of the DS 



is PddW ~ *( JJ l)/ 8rw 2i)Pi2 W + c - c - Therefore the 
contribution of the fast-oscillating terms to the observed 
probability R(t) is < (JJ 2 ) /2^ 1 )Repi 2 (0). It is small 
compared to the measurement error ~ J 2 /wfi that we 
want to overcome. 



({), 



APPENDIX B: DYNAMICS OF THE 
TWO-EXCITATION STATE 

In the presence of two excitations, along with the two- 
particle matrix elements p a p{t) = (ag(£)a a (£)) the dy- 
namics of the system is characterized by the four-particle 
matrix elements p a pys(t) = (aJ(£)a^,(£)a J g(£)a Q (£)). For a 
two-qubit system coupled to a DS each Greek subscript 
runs through 3 values, and the dynamics is described by a 
set of 18 linear differential equations that follow from the 
equation of motion for the density matrix in the operator 
form ((4]): 9 equations for p a /j and 9 equations for Pa/3-yS- 
The number of independent matrix elements follows from 
the permutation symmetry S ^ 7, f3 ^ a (a permutation 
is accompanied by a sign change) and the commutation 
relations for the fermion operators a a , a^ a . The system of 
equations for the four-particle matrix elements is closed: 
PapyS are expressed in terms of each other, whereas their 
evolution affects the two-particle matrix elements p a p. 

We will again use the adiabatic approximation to study 
the dynamics and will start with the matrix elements 
Pap-yS- Formally, their time evolution is determined by 
the eigenvalues of the system of equations for p. Four 
of these eigenvalues correspond to fast oscillations with 
imaginary part ss ±u>2i- In the limit J — > the corre- 
sponding eigenvectors Papyg have one of the subscripts 
equal to 2, whereas other three subscripts are two 1 and 
one D or two D and one 1. Among the eigenvalues that 
describe "slow" motion there are two w — V and two 
w — 2r. In the limit J — > the corresponding eigen- 
vectors arc pD2i2,Pi2D2 and p rm DD {n = 1,2). Only one 
eigenvector controls the evolution of the four-particle ma- 
trix elements on times longer than the leading term 
in this eigenvector is pi22i- 

Over time t > T^ 1 the matrix elements p a /3-yS reach 
their quasi-stationary values which are determined by 
Pi22i(£)- They can be found by disregarding pap^s in all 
equations except for the equation for pi22i. This gives, 
in particular, 



AD22l(£) 



iJp 
~2~F 



P1221W, P1221 » -W1P1221, (Bl) 



whereas p\D2\ ~ i(J Jd/ ^>2i^)pi22i- It is seen from 
Eq. (|Blj) that all slowly varying components of p a /3~f5 de- 
cay to zero over time ~ W± . This makes sense, since 
this is the time over which the occupation of state {tpn^) 
decays, after which there remains only one excitation in 
the system. The fast-oscillating components of Pa/3-yS de- 
cay even faster, over time 

The signal R(t), Eq. ([5]), is determined by the two- 
particle matrix element pdd (£) ■ One can show that, 
in the time range £ > T" 1 , the major result of the 
interaction between the excitation on the two-particle 
matrix elements is that pD\{t) is incremented by sa 
(iJA/T)p D2 2i(t) and the slow part of PD2(t) is incre- 
mented by w — (JA/w 2 i)/Oiz52i(£). The first of these cor- 
rections drops out of the expression for pdd, whereas 
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the contribution from the second is small compared to 
(J/w 2 i) 2 - Therefore, to the accuracy we are interested 
in, the interaction between excitations does not affect 
measurements in the proposed approach. 
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